WWM]SSf:1flNSRIRflM Fax: 7132668510 



Dec 6 '04 11:05 P. 02 



RECEIVED 
CENTRAL FAX CENIER 

DEC 0 6 20M 



ore 15125 

Real-time Estimation of Resistivity Anisotropy Using Array Lateral and Induction Logs 
Michael A. Fr«nkel and ingo M. Geldmacher/ Baker Hughes 




7KI p^pv «»■ pnpaivd pmtfilaftm St Ihi 2lim QiM»m 
HOitflon. TwM, U AA. May 2flOa. 



eoRtMM n an Msna wtoAM ey Ita ausioi^s). CoAtwt* or ttw ps^. sa 

' ~ " * ^^taJbi«t» 



iiiMnfj. Nw* eia »ae«t rvtimM ty Oifchow Tthreiwy Omftetnco and 
evi«0bR by Aft «U3)0rt«^ malBrM, at mwrtcd. dM& ml nM««*«i«y 
aaOlon «« B" Ofhtio™ T^enwogy oonterenw or orcaM. aectionkc n 
tffjTa..<t,y.^ ^ -;>yafln flf any part ffih taflf tof oorwnttfcid Ourpo*W ^waiout tt» wnBcn 

n Btntmet «l n«t more Bon 300 ovstnnlons may not b* eopi^- Tha 



Abstract 

Conventional resistivity tools often miss hydrocarbon pay 
zones in thinly laminated sand-shale sequences. iDcorpcxation 
of fonmtioo resistivity ai^isotropy into & petrophysical model 
allows ior a significantly more accvmte means to cstioiate 
hydrocarboa reserves at low-resistivity reservoirs. To 
correctly detenmne tide anisotropy distribution srornid Ibe 
borehole, an intetpreter must use flie mutti-conxponent 
inductiDn logs acquired with the appiicaKon of the new 
3D£X tool. 

If the imilti-component inducfion measurements are not 
mailable, an fnteipreter should consider some alternative ways 
to determine aoisonopic propoties of the fomxation. 

In this paper, we propose two new methods to estimate 
resistivity anisotropy using vast legging data from the Array 
Lateral Log (High-Definition Lateral tog ■ HDLL) and Array 
Induction Log (High-DeEciition Induction Log -HDIL) tools: 



■ The first method is based ODi a sequential 
interpretation of the HDO. and HDLL data. It does 
not require layer selection ftir fonnation model 
definition. It is very fast because it uses sequential 
induction aiod g^^^^^ processing as wdl as 
short iiwez^ioo windows. 

■ The second method utilizes boreholc-coirected lateral 
resistivity (LR) logs and vertical rcsohition matched 
(VRM) focused logs. The method does not require 
application of inversion-based processing and 
provides anisotropy values at every logging point. 
Tliis approach is extresnely fest. 

Wc evaluate and compare the anisotropy iDtexpretailon lesolts 
derived from the application of these methods. We \ise a 
150-foot portion of a single data set acquired in a vertical 
OfBhore well in the Me^teiranean region. This unique 
logging data set contains various different anay lesisttvity 



logs. Wc show that both methods supply an interpreter with 
anisotropy cscitnates whhin time limits in^TOSed by real-time 
well*site processing. 

Background 

A significant part of the world's estimated hydrocarbon 
reserves is contained in thiniy-bedded, low-resiStivity 
formations containing hydit>€arboQ producing laminar sands. 
In these fonnations, the transvetse (vertical) resistivity, Rv» 
perpendicular to the bedding plane is greater than the 
longitudinal (horizontal) resistivity^ ^ parallel to the bedding 
plane. This dependence of resistivity on the direction of 
cuirent flow is called anisotropy. 

The anisolTOpy coefiScient is dcfmed as a ratio Of K, 
over K^; 

3 = ^* 

The xradiUonal laminar shaiy sand saturation equations are 
based on a simple parallel condactivity model that Qipieally 
results in significant underestimation of - iQ^droeaibonft- 
in-plaoc. The resistiwty anisotn^ based Inteipretation 
rocthodsj when combined whh peiropfaysical analysis, lead an 
interpreter to a better identification and a noore accurate 
quantillcaxton of hydioeaxbon reserves.^ 

An Bocuratt determinafi«i of f<»inatioQ anisotropy could 
be achieved via (he applicatioso of unique miilti-eomponcnt 
induction (3DEX) technology. This approach was developed 
to provide both improved identification and quantificatioa of 
hydrocarbon*. Our 3DEX tool provides all necessary data to 
accurately compute both horizontal (R<0 and vertical (R^) 
resistivities.''* Bohanced shaly sand reservoir characterization 
models utiHzlng multi-component induction data lead 
interpieieia to a significant increase ia the estimates of 
hydrocaibons^in-place. * 

If the multi-component induction measurements are not 
available, or if use of this instrument under given 
environmental conditions is not recoimnendable^ an interpreter 
should consider alternative ways to determine the anisotropic 
properties of the fonnatxon. 

The effect of foirnation atiisotropy on conventionai 
resistivity measurements have been well studied in numerous 
publications.*''' In near-vertical wells, logging tools such as 
Dual taterolog (PIX) or Airay Induction general^ currents 
primarily paraUel to the fonnation bedding layers, and are thus 
itot senative to changes of the formation resistivity in die 
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vertical diitction. Therefore, even a jgint interpretation based 
on DLL-typc, galvanic logs and array induction-type logs** 
cannot gencrtJIy provide reliable anisotropy estimates. 

On the other hand, the lateral-type, unfocused log contains 
information related to both horizontal (Rj,) and vertical (R,) 
resistivities. Moreover, If wc neglect the effects of both the 
borehole and the invaded zones, an anisotropic bed can be 
replaced by fin equivalent isotropic bed describable with the 
following parameters^: 

h^- VI h (1) 

In (1) and (2), h and Rh are the tnic fonnaelon thickness 
and resistivity, respectively; \^ and are equivalent 
fonnation thickness and resistivity, respectively; and X is the 
Bnisonropy ratio. 

This besucifUl equivalence imtDediatcly suggests several 
alternative q[iproacbes to estimate fonnation aAlsotropy. 
These inclnde using an unfocused lateral log in combination 
with additional logs providing qccuhUjo infofmation about true 
formation thic^ess and horizontal rcdsttvity. 

One approach is based on the joint 2-D inversion of lateral 
and indwUoa logs.* Tor the first step, the author suggests 
recovering the true layer thickness h and its horizontal 
rtSistivity by using an induction log. For dae second Step, 
and to calculate the anisctropy ratio, the author uses tiie results 
of the first step combined witibi results from a lateral 
log invcxTion. 

The method* was tested on three-layer models using only 
conventienAl latwal and induction synthetic logs. There are 
several shortcomings to ibis approach. Firstly, the practicality 
of this method will require careful layer picking. Secondly, to 
detemiine an accurate Bh, an interpreter must cany out the 
inversion of an induction log. Thirdly, the limited number of 
eoflventional galvanic and induction logs available when using 
this method may not provide enough data for a more accurate 
estimation of fiumation anisotropy. 

In the following sections of this paper, we first briefly 
describe the fiimily of available array resistivity tools. Next, 
wc describe our two new methods for calculating formation 
fcsbtivitjr anisotropy. Finally, we test these iaethods on a 
single weU logging data set condscing of several different 
aiiay logs. We then dcmomstrste thai these methods can 
provide real-time anisoinw estimates even in difficult 
field situations. 

Array Resistivity Tools (HDIL, HDLL, 3DEX) 
The HDIL h an array induction tool that collect data at 
multiple ftcquencies and vaiioms tmsniitter-rcceiver 
spacings.' This tool utilizes only vertical tjansmitter and 
receiver anay coils (the coil centers arc m line with the tool 
axis). Therefore, the induced current in near-vertical wells is 
dominated by io horizontal component, and the induction data 
contains infonnatlon related to alone. 

Existing fbeusing algorithms arc designed to convert the 
HDIL measurements into ao-called Vextieal Resolution 
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Matched (VKM) logs. In shaly-sand lyne fcscrvoirs, the 
deep-reading VRM curves read close to Rh. 

The HDLL is an array-type, unfocused galvanic logging 
tool. It has a single current injection electrode. The tool 
mcastires an array of the absolute potentials and the fir^it 
differences at selected rectnver electrodes. On the basis of (he 
measured finx dtfibrcnccs, a set of lateral Fe«istivtty (RL) logs 
are compuced. 

In near-vcrtical wells drilled with conductive drilling 
muds, the HDLL instrument provides vertical resohition and 
resistivi^ invasion profiles simitar to ihoso provided by array 
induction tools in oil-baaed mud." Since the injected cuirent 
has both horizontal and vertical components, the acquired data 
contains information related to both horizontal 
vertical (R») resistivities. 

The 3DEX tool employs sets of Z direction coils coaxial 
witb the instrument id addidon to orthogonally mounted X and 
Y direction coil arrays. This setting provides measurements of 
five components: Hn. THyp H^y, and Km ^ 

seversal frequencies- 
Inversion processing of 3DHX induction data then allows 
the computation of both horizontal (R^^ and vertical (Rv> 
resistivities. In mm, this icads to an accurate determination of 
the fbrmation resistivity amsotropy ratio, X. The incorporation 
of such induction data intexpretation results into an enhanped 
shaly-sand, tensor resistivity petrophysical analysis leads to 
much reduced evaluation uncertainties- The inclusion of 
multi-componeiit induction data may thus result in a 
significant increase in calculated l^drocaibon^in-place 
Tcseives over estimates obtained witibi convendonal 
methodologies.' 

Anisotropy Estimatioi^ by Inversion of the Array 
Lateral and Induction Logs 

It is well-known that the deep-reading, borebolc'COtrectcd 
VRM curves read R|| quite accurately in shaly sands."* 
Thcmfore we do not have to run the array induction data 
inversion but assumo that R|, = VRM. To make the anisotropy 
estimation process automatic and robust; we suggest squaring 
the devest reading VICM curve based on a ccmstont 
layer thickness. . . 

We tibien use this model as a starting point and proceed to 
invert the sdected lateral logs. This 2-D invenion allows ^X3. to 
correct tiie data for borehole, invasion, and shoulder bed 
effects. As ou^ut of this process, we obtain parometers of 
equivalent isotropic layers. 

Finally, we can detemaine X from equation (2) by usbg Rn 
derived firom the azray induction loga, and R^, derived fiom 
the array lateral logs- 
It should be noted that application of a short ovcriapping 
inversion windows rnakes the i>roce$3 extremely fast **>an 
one ininute/lOa feet of data on SunBladclOO). 

Anisotropy Estimation with Correctad Array Lateral 
and Induction Logs 

The method described in the previous section can be 
simplified and accelerated. As before, we get Ri from the 
deep-reading, borehole-corrected VRM curve. Then, instead 
of tunning the 2-D inversion of the lateral logs, we perform 
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cbart-4>9sed log correction for the borehole and, if necessary, 
invasion effects. 

To esttmale (he formation anisotropy, we use the following 
sppioximate fonnula: 

The calculations arc perfonned at each logging point. This 
provides % oonttmious smiwhopy output log. The method 
should thus fhmish reasonable anisotropy esdmatcs in thick 
sand/shale fannadons where Che borohole-conectcd RL logs 
are not &igtiJflcamly influenced by sigmfieant shoulder 
bed effects. 

In the next section we present a case study that illustrates 
Oie practiGa] application of out proposed new methods. 

A Case Stuffy 

MuW-compOnent induction (3DEX) and convenrional wireline 
array lesbtivity logs (HDIL and HDLL) were acquired to a 
. vertical oflfthore well in the Mediterranean region. The well 
-was drilled xatct a deep xnaxine ftUblditc scquimce. It 
encountered significant volumes of dunly bedded, laminar 
silty shales as well as high porosity gas sands deposited over 
very high quality, massive channel and (uHridite ^ sands. 

In Fig. 1, we present the 2-D inveraion tcsults for a 
150-fooc interval. 600 -750 fl (depths arc relative). The left 
track (1) shows the Caliper and GR logs. 

Traclcs (2-3) show Rh and R«, cstimated ftom HDEL and 
HDLL. ntspectively. using the two methods presented in this 
p^er. The separation between Rd and R«, is evident in both 
right tracks, and indicates the degree of fonnation resistivity 
a^sotropy within this interval. 

It should be noted that inversion-based processing with 
short, sliding down windows takes less than 1.5 miixute on a 
SuriBIadclOO computer. 

In Hg, 2, we present the results of anisotropy esticfialion 
with our methods. We also compare ilicsc results with the 
tBftrcnw anisotropy model wc derived with the accurate 
3D£X/iiE>IL data invexslon method.* The step-wise curves 
wcfo smoothed out and we pieseat tbem as coodnuous logs. 

The left two tracks show die Caliper, GR, Dcnshy, and 
Neutron logs. Track (3) shows a comparison between the 
anisotropy factor, AFl, derived with ow first method versus 
tfac i^fereaco anisotropy facnar, AF2, derived from application 
of the mnlti-componeni induction data. Wc recall that the firtt 
method is based Ott HDlL-focuaed data ami HDLL 
2-D invetsioEL 

Track (4) shows the comparison between the anisotropy 
fector. AFl, versus the Anisotropy factor, AF3, derived from 
our method based on the correction of the HDIL-focuscd logs 
and HDLL lateral logs. 

These comparisons show a reasonable match between the 
anisotropy fectors, AFl and AF2, along tha entire interval, 
except in the 710-725 feet section, where the reference 
anisotrxipy gives higher values. The match between the 
anisotropy factons AFl and AF3 is worse than between 
anisotropy factora AFl and AF2. This could be attributed to 



the influence of the uncorrected shoulder bed effect on the RL 
log, which might well affect the AF3 estimates. 

. Our results provide us with a practical validation of die 
pr^osed mcdiods, and lead to the following ooncluaioias. 

Conclusions 

Wc find the following: 

■ The multi-ooirtponent induction (3DE70 technology 
provides the most accurate method to determine 
formation resistivity anisotropy. 

• Joint interpretation of the Array Induction (HDIL) 
and Array Lateral Log (HDLL) data can be used as 
aa alramative way to estnnate formation 
resistivity anisotropy. 

• The anisotropy interpretations obtained with our 
proposed me^ds are in reasonable agreement with 
an independent anisotropy in terpr etation based ort the 
direct inversion of multi-compoorat Indocdon data 
acquired in die same welL 

• The methods are fully automatic and can be used in 
situations when superior multi*componcnt induction 
data is not available. 

■ The methods could be used fbr the lepiocessing of 
previoosly acquired tnduciioa and lottitd logs. 

Future developments should include: 

- Expansion of the methods to mote complicated, 
invaded-type formation models. 

• More rigorous evaluation of the methods for both 
syxithetic models and «»r addiidonat field cases. 
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Nomenclature 

k - True thicJcncss of anisotropic l<tyer, m 

h^^ ■= Thickness of equivalent isotropic layer, m 

Rf, — Horizontal resisiivi^ of anisotropic fomtation, OhftMn 

= Resistivity of equivalent isotropic formation, Ohm-m 
R^ = Vertical resistivity of anisotropic formation, Ohm-m 
RL = Laterat Rexisttvfty curve, Ohm-m 
VRM = Vertical Resolution Matched curve, Ohm-m 
X = Anisotropy ratio 
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